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ABSTRACT 

Aims. The aim of this study is to understand the chemical conditions of ices around embedded young stellar objects (YSOs) in the 
metal-poor Large Magellanic Cloud (LMC). 

Methods. We performed near-infrared (2.5-5/um) spectroscopic observations toward 12 massive embedded YSOs and their candidates 
in the LMC using the Infrared Camera (IRC) onboard AKARI. We estimated the column densities of the H 2 0, CO2, and CO ices 
based on their 3.05, 4.27, and 4.67^m absorption features, and we investigated the correlation between ice abundances and physical 
properties of YSOs. 

Results. The ice absorption features of H 2 0, CO2, I3 C02, CO, CH3OH, and possibly XCN are detected in the spectra. In addition, 
hydrogen recombination lines and PAH emission bands are detected toward the majority of the targets. The derived typical CO2/H2O 
ice ratio of our samples (~0.36 ± 0.09) is greater than that of Galactic massive YSOs (~0.17 ± 0.03), while the CO/H 2 ice ratio 
is comparable. It is shown that the CO2 ice abundance does not correlate with the observed characteristics of YSOs: the strength of 
hydrogen recombination line and the total luminosity. Likewise, clear correlation is not seen between the CO ice abundance and YSO 
characteristics, but it is suggested that the CO ice abundance of luminous samples is significantly lower than in other samples. 
Conclusions. The systematic difference in the C02 ice abundance around the LMC's massive YSOs, which was suggested by previous 
studies, is confirmed with the new near-IR data. We suggest that the strong ultraviolet radiation field and/or the high dust temperature 
in the LMC are responsible for the observed high abundance of the C0 2 ice. It is suggested that the internal stellar radiation does not 
play an important role in the evolution of the CO2 ice around a massive YSO, while more volatile molecules like CO are susceptible 
to the effect of the stellar radiation. 
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1. Introduction 



The infrared spectra of embedded young stellar objects (YSOs) 
show absorption features originating in various ice species (solid 
state molecules, e.g., Gerakines et al. 1999; Gibb et al. 2004; 
Boogert et al. 2008). It is believed that a large amount of heavy 
elements and complex molecules of the interstellar medium 
(ISM) are preserved as ices in the dense and cold regions (n# > 
10 4 cm 4 , T ~ 10 - 20K), such as an envelope of a deeply embed- 
ded YSO (van Dishoeck & Blake 1998; Boogert & Ehrenfreund 
2004). Since the star formation proceeds in the region of dense 
molecular clouds, these ices dominate the chemical evolution of 
YSOs. Absorption profiles of ices are known to be sensitive to a 
chemical composition and a temperature of dust grains, and ices 
are important tracers of the thermal history of circumstellar envi- 
ronments of YSOs (e.g., Pontoppidan et al. 2008; Zasowski et al. 
2009). Thus, investigating the compositions of ices as a function 
of physical environments is crucial for understanding the chem- 
ical evolution of YSOs and is one of the key topics in the cur- 
rent astrochemistry. Ices are also detected toward solar system 
objects such as comets, icy satellites, and Mars. Interstellar ices 
are thought to be taken into planets and comets as a result of sub- 
sequent planetary formation activities (Ehrenfreund & Schutte 
2000). Therefore the chemical evolution of ices around YSOs is 



also of interest in terms of chemical conditions of the planetary 
formation. 

Observations of ices around extragalactic YSOs are one of 
the challenges in current ice studies (Shimonishi et al. 2008; 
Oliveira et al. 2009; Seale et al. 2009; van Loon et al. 2005, 
2010a). So far, infrared spectroscopic observations of extra- 
galactic embedded YSOs are still limited by observational dif- 
ficulties, and their circumstellar chemistry is poorly understood. 
However, it is probable that different galactic environments (e.g., 
metallicity or radiation field) could affect the properties of cir- 
cumstellar material. Thus detailed investigations of extragalac- 
tic YSOs provide us important information for understanding 
the diversity of chemical conditions of YSOs in the universe. 
Thanks to the progress of infrared space telescopes, such as 
AKARI (Murakami et al. 2007) and Spitzer (Werner et al. 2004), 
we are now able to extend the study of ices around YSOs to ex- 
tragalactic objects. 

The Large Magellanic Cloud (LMC), the nearest irregular 
galaxy to our Galaxy, offers an ideal environment for studying 
extragalactic star formations since individual YSOs can be iden- 
tified with reasonable spatial resolution (l"~0.25pc) thanks to 
its proximity (~50kpc; Alves 2004). A unique metal-poor envi- 
ronment of the LMC (Z ~ 0.3Z o , Luck et al. 1998) results in 
a generally high UV radiation field over the galaxy, and these 
factors can affect the chemical conditions of circumstellar ma- 
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terials. The nearly face-on viewing angle of this galaxy (~35 
degrees, van der Marel & Cioni 2001) allows correlation studies 
of individual YSOs and the ISM. Because of these advantages, 
various types of surveys have been performed toward the LMC, 
which provide information of both ISM and stars over a wide 
wavelength range (e.g., Zaritsky et al. 2004; Meixneret al. 2006; 
Kato et al. 2007). 

Shimonishi et al. (2008) spectroscopically confirmed seven 
massive YSOs in the LMC based on the AKARI LMC spectro- 
scopic survey (Ita et al. 2008) and reported the detection of the 
H2O and CO2 ices toward these objects. Oliveira et al. (2009) 
analyzed the 15.2yum feature of the CO2 ice for 15 embedded 
YSOs in the LMC and discuss the chemical properties of the ice 
mantles. The H2O and CO2 ices are ubiquitous and are the major 
components of interstellar ices. Based on the low-resolution (R 
~ 20) spectra obtained in Shimonishi et al. (2008), we derived 
a typical CO2/H2O ice ratio and showed that the ratio around 
the LMC's massive YSOs (0.45±0.17) is higher than previously 
reported toward Galactic massive YSOs (0.17±0.03 Gerakines 
et al. 1999). The result suggests the chemically different nature 
of YSOs in the metal-poor galaxy. 

An accurate determination of column densities of major ice 
species is crucial when discussing the variation of ice abun- 
dances between the objects. However, the uncertainties in the de- 
rived column densities from our previous low-resolution spectra 
are large, and the number of samples for which the column den- 
sities of major ice species are determined is still small. In addi- 
tion, the features of relatively minor ice species, such as CH3OH 
and CO are difficult to identify by the low-resolution spectra. 
Investigation of these minor ices should help in constraining the 
processing of ices in a YSO envelope (e.g., Dartois et al. 1999; 
Thi et al. 2006). Thus follow-up observations with higher spec- 
tral resolution are required to further understand the chemical 
conditions of extragalactic YSOs. 

In the present study, we used the higher spectral resolution 
(R ~ 80) mode of the Infrared Camera (IRC, Onaka et al. 2007) 
onboard AKARI to observe massive embedded YSOs and their 
candidates in the LMC. The correlation between the chemical 
properties of ices and the YSO characteristics, which could not 
be discussed in Shimonishi et al. (2008) due to the large uncer- 
tainties of the observed spectra, is discussed in this paper. In ad- 
dition, the improved calculation method of the ice column den- 
sity by using the curve-of-growth is presented in an Appendix. 

2. Targets 

The sample of massive YSOs presented in this study is mainly 
selected from Shimonishi et al. (2008), which spectroscopi- 
cally discovered several massive YSOs in the LMC by us- 
ing a combination of the low-resolution, near-infrared spectro- 
scopic survey and the color-color diagram. The H2O and CO2 
ice features toward these targets have already been reported in 
Shimonishi et al. (2008), and this indicates the embedded na- 
ture of these YSOs. In addition, we include several YSO "can- 
didates". These objects were selected by the same photometric 
criteria described in Shimonishi et al. (2008) and their infrared 
spectral energy distributions (SEDs) are red enough to be iden- 
tified as a YSO. But absorption features of ices seen in their 
low-resolution 2-5yum spectra are too weak to definitely clas- 
sify them as YSOs, so they are not included in the YSO list 
of Shimonishi et al. (2008). Some other targets were selected 
based on the YSO candidates catalog presented in Whitney et al. 
(2008), which performed photometric selection of YSOs using 
the database of Spitzer SAGE project (Meixner et al. 2006). We 



selected relatively bright and red objects from their catalog. Van 
Loon et al. (2005) observed IRAS05328-6827 with the Infrared 
Spectrograph (IRS, Houck et al. 2004) onboard Spitzer and the 
Infrared Spectrometer And Array Camera (ISAAC) at the Very 
Large Telescope (VLT), and report the presence of ice features 
in its near- to mid-infrared spectrum. This object is also added 
to the present targets. 

To understand the nature of the targets, we carried out the 
SED fitting of each YSO by using the Online SED Fitter 1 
(Robitaille et al. 2007). The photometric dataset used for the 
SED fitting are obtained from the available database and recent 
YSO catalogs (Zaritsky et al. 2004; Meixner et al. 2006; Kato 
et al. 2007; Ita et al. 2008; Gruendl & Chu 2009), and 2-24//m 
data were typically used as the input data. The estimated total 
luminosities and masses of our targets are ~4-34OxlO 3 L and 
12-47M Q , respectively. Thus the targets in the present study are 
massive YSOs. Details of all the targets are listed in Table 1. 

3. Observation and data reduction 

Spectra of 2-5/im presented in this paper were obtained 
as a part of AKARI post-helium open time program "Ices 
Around Extragalactic Young Stellar Objects" (IEYSO, PI: T. 
Shimonishi) and AKARI Director's Time (DT) observations. 
AKARI is the first Japanese satellite dedicated to infrared as- 
tronomy launched in 2006. The telescope system and the scien- 
tific instruments onboard were cooled down by the liquid helium 
with mechanical coolers. The AKARI cryogen boiled off on 2007 
August 26, 550 days after launch. All of the data presented here 
were obtained after the exhaustion of liquid helium. 

The slit-less NG spectroscopy mode with the Np point 
source aperture was used to obtain spectra between 2.5 and 5/im 
with a spectral resolution of R~80. The wavelength accuracy of 
the spectroscopy is estimated to be ~0.01-0.015//m. The spatial 
resolution of a point source is approximately 5-8 arcseconds. 
Overlapping of dispersed spectra with other objects is a serious 
problem for slit-less spectroscopies, which makes extraction of 
reliable spectra difficult. The use of the Np point source aperture 
reduces the effect of spectral overlapping and enables spectra to 
be extracted from the objects located in crowded regions (see 
Ohyama et al. 2007, and AKARI Data Users Manual for details 
of the IRC spectroscopy). 

The basic spectral analysis was performed using the standard 
IDL pipeline prepared for reducing AKARI IRC post-helium 
mission data. Raw data were converted to dark-subtracted and 
linearity-corrected frames in the pipeline. We used a newly cal- 
ibrated spectral response curve to analyze the present spectra. 
At that time, appropriate flat field data were not available for 
post-helium Np data. However, because of the faintness of the 
natural background and the small pixel-to-pixel variation of the 
flux at the Np aperture, there are no significant differences be- 
tween the spectra with and without flat-fielding (AKARI Data 
Users Manual ver. 1.3). In fact, spectra processed without the 
flat-fielding is better in terms of S/N because of the relatively 
low S/N of the flat field data. Therefore, we did not perform the 
flat-fielding in the pipeline. Four to seven pixels were integrated 
in the spatial direction to extract the point source spectrum. Each 
pixel was smoothed in the dispersion direction across the width 
of three pixels to suppress the effect of bad pixels. Background 
signals of each source were estimated from the signal counts of 
the adjacent region of the target spectrum and subtracted. 

1 The Online SED Fitter is available at 
http://caravan.astro.wisc.edu/protostars/sedfitter.php 



Takashi Shimonishi et al.: Ices around Embedded YSOs in the LMC 



3 



Table 1. Targets of the present observations 



Number 


ID 


Obs. ID 


Obs. Date 


RA[J2000] 


DEC[I2000] 


Reference 


ST1 


IRAS05400-7013 


5126560 


2008 Apr 12 


05:39:31.15 


-70:12:16.8 


1,2 


ST2 


NGC 1936 


5126562 


2008 May 25 


05:22:12.56 


-67:58:32.2 


1 


ST3 


a 


3410003 


2008 Dec 11 


05:25:46.69 


-66:14:11.3 


1, 3 


ST4 


IRASF05148-6715 


3410004 


2008 Dec 27 


05:14:49.41 


-67:12:21.5 


1,3 


ST5 


IRAS 05311-6836 


5126561 


2008 May 1 1 


05:30:54.27 


-68:34:28.2 


1, 2 


ST6 


05394112-6929166* 


5126563 


2008 Apr 19 


05:39:41.08 


-69:29:16.8 


1 


ST7 


IRAS05240-6809 


5126567 


2008 May 23 


05:23:51.15 


-68:07:12.2 


1,3 


ST8 


MSX LMC 1786 


5126573 


2008 May 2 


05:37:28.17 


-69:08:47.0 


3,4,5,6 


ST9 


IRAS05237-6755 


5126568 


2008 May 26 


05:23:35.53 


-67:52:35.5 


2,4 


ST10 


MSX LMC 1229 


5200384 


2008 Dec 21 


04:56:40.80 


-66:32:30.4 


2, 4, 5 


ST11 


IRAS05270-6851 


5126571 


2008 May 1 1 


05:26:46.63 


-68:48:47.1 


2,5 


ST12 


IRAS05328-6827 


5126566 


2008 May 12 


05:32:38.59 


-68:25:22.2 


3,6 



"The source is in a cluster, *2MASS ID 

References. (1) Shimonishi et al. (2008), (2) Seale et al. (2009), (3) Oliveira et al. (2009), (4) Whitney et al. (2008), (5) Gruendl & Chu (2009), 
(6) van Loon et al. (2010a) 



Table 2. Details of detected ice features 



Molecule 


Wavelength 


Vibration Mode 


A" 


Reference 




[pan] 




[10~' 7 cm molecule -1 ] 




H 2 


3.05 


O-H stretch 


20 


1 


CILOH 


3.53 


C-H stretch 


0.76 


2 


co 2 


4.27 


C-O stretch 


7.6 


1 


13 C0 2 


4.38 


13 C-0 stretch 


7.8 


1 


XCN 


4.62 


CN stretch 




3 


CO 


4.67 


C-O stretch 


1.1 


1 



"Band strengths of absorption features 

References. (1) Gerakines et al. (1995), (2) D'Hendecourt & Allamandola (1986), (3) Schutte & Greenberg (1997) 



4. Results 

4.1. Observed spectra 

Observed spectra for all the objects are shown in Fig. 1 and de- 
tails of the detected ice features are summarized in Tables 2 and 
3. The absorption feature of the 3.05/im H2O ice is detected to- 
ward all the objects except ST6 and ST9. The profile of the H2O 
ice feature is resolved, and this permits more accurate determina- 
tion of the column density than in Shimonishi et al. (2008). One 
object (ST8) shows a structured profile of the H2O ice, which is 
often linked to the thermal processing of the ice mantles (Oberg 
et al. 2007). The water-ice band of ST6 is seen in the spectrum of 
Shimonishi et al. (2008), but it is not detected in the present spec- 
tra owing to the slightly poor sensitivity of the present observa- 
tion mode. The 3fim region of ST9 overlaps prominent emission 
components, which makes it difficult to identify the underlying 
absorption bands around 3/im. A broad depression is seen in the 
4.3 - 4.7yum region of ST7. It can probably be attributed to a 
combination mode of the H2O ice, which typically peaks around 
4.5//m (Gibb et al. 2000). 

The absorption feature of the 4.27//m CO2 ice is detected to- 
ward all the objects. Some objects show the absorption features 
of 3.53//m CH3OH ice, 4.38/mi 13 C0 2 ice, 4.62^111 XCN, and 
the 4.67//m CO ice. Since the CO2 and CO ice features are not 
completely resolved in our spectra, the widths of the observed 
features are larger than those of their typical laboratory spectra. 
The detections of the 13 C02 ice and XCN features are still tenta- 
tive because of the unresolved absorption bands of gaseous CO, 
which sometimes appear in the wavelength region of these fea- 
tures (Chiar et al. 1998; Boogert et al. 2000). The non-detection 



Table 3. Detection of ice features in our YSO spectra 



Number 


H 2 


co 2 


CO 


CH30H 


IJ co 2 


XCN 


ST1 


V 


V 


V 








ST2 


V 


V 




V 


? 




ST3 


V 


V 


V 








ST4 


V 


V 


V 








ST5 


V 


V 


V 


V 






ST6 




V 


V 






? 


ST7 


V 


V 


V 




V 


? 


ST8 


V 


V 


V 




V 




ST9 




V 


V 




? 




ST10 


V 


V 


V 






? 


ST11 


V 


V 




? 






ST12 


V 


V 




? 







Note. V = c l ear detection, ? = tentative detection, — = no detection. 



of the 13 C0 2 ice in ST10, on the other hand, gives a rough lower 
limit of 13 C02/ 12 C02 of about 95. This is in the same range as 
the isotopic ratio for Galactic YSOs (Boogert et al. 2000). These 
minor features are labeled in Fig. 2 and summarized in Table 3. 

Several emission features that were not resolved in the low- 
resolution spectra of Shimonishi et al. (2008) are detected in the 
present spectra. They are hydrogen recombination lines (2.62//m 
Br/3, 3 .05//m Pf e, 3 .29yum Pf6, 3.74/um Pfy, 4.05//m Bra, 4.64//m 
Pf/T) and polycyclic aromatic hydrocarbon(PAH) emission bands 
(3. 3/im and 3.4/im). Detection of these emission features allows 
detailed discussions about the correlation between ice abun- 
dances and YSO characteristics (see §5.2). 
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4.2. Spectral fitting and calculation of ice column densities 

To subtract the continuum, we fitted a polynomial of the second 
to fourth order to the continuum regions and derived the optical 
depth spectra. The wavelength regions used for the continuum 
were set to avoid a region with emission or absorption features, 
and they were typically 2.5-2.7 /mi, 4. 1-4.15 fim, and 4.85^4.95 
/zm. The absorption features of ices except H2O were not com- 
pletely resolved with our spectral resolution. Therefore, to cal- 
culate ice column densities, we used different methods for the 
H2O ice and the other ices. 



4.2.1. H 2 ice 

The spectral resolution of the present spectra is sufficient to re- 
solve a broad feature of the H2O ice, and the true optical depth 
can be measured. Thus a direct comparison of laboratory spectra 
is possible. We fitted laboratory ice profiles to the observed spec- 



tra by ax 2 minimization method and derived the column density 
from the equation 

N= J rdv/A, (1) 

where A is the band strength of each ice feature (in units of cm 
molecule -1 ) as measured in the laboratory, N the column density 
in crrT 2 , v the wavenumber in units of cirT 1 , and r the optical 
depth. The integration is performed over the H2O ice absorption 
band. The wavelength range used for the fit is set to avoid the 
3.3yum PAH emission band and the bottom of the H2O absorp- 
tion feature, which typically suffers from the weak background 
emission. Typically, 2.7-2.9/im and 3.1-3.2/im are used for the 
fitting. A single laboratory ice profile of the pure H2O ice at 10K, 
which is expected to exist in the cold envelope of YSOs, is fit- 
ted to the observed spectra The laboratory spectra of the H2O 
ice was taken from the Leiden Molecular Astrophysics database 
2 (Ehrenfreund et al. 1996). We adopted the band strengths of 

2 data available at http://www.strw.leidenuniv.nl/ lab 
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Fig. 2. Observed spectra of sample YSOs after the continuumd subtraction. A vertical axis represents the optical depth. The positions of the 
detected features are labeled. (Left) 2.6 - 3.8/um region. A solid line represents the fitted laboratory spectrum of the H 2 ice (pure H 2 at 10K, 
Ehrenfreund et al. 1996). 2.9 - 3.1yum (bottom of the H^O ice feature) and 3.2 - 3.4yum are removed from the fitting wavelength. (Center) 3.9 - 
4.4yum region. A solid line represents a fitted Gaussian profile for a measurement of the equivalent width of the CO2 ice feature. (Right) 4.4 - 
4.8yiim region. A dashed line indicates unresolved CO gas bands and the possible XCN feature, which were subtracted when fitting the CO ice 
feature (solid line). 



the H2O ice band as 2.0xl(T 16 cm molecule 1 (Gerakines et al. 
1995). The derived column densities are listed in Table 4, and 
results of the fitting are shown in Fig. 2. 

4.2.2. C0 2 ice and CO ice 

In case the spectral resolution is insufficient to resolve the ab- 
sorption features, we cannot derive the "true" optical depth from 
the spectra and we can not use Eq. (1) to calculate column den- 
sities. Instead, we use the equivalent width of the absorption fea- 
tures to calculate column densities since it does not depend on 
the spectral resolution. The equivalent width (EW) is defined as 

EW = f F norm dA, (2) 



and F norm is the normalized flux, F n i, s and F c are the observed 
and continuum flux, and A is the wavelength. A Gaussian profile 
with the fixed central wavelength was fitted to the CO2 and CO 
absorption features to measure the equivalent width. To convert 
the equivalent width to the column density, we need to calcu- 
late the "curve-of-growth" of each absorption feature. The de- 
tailed discussion about our method is described in Appendix. A. 
To calculate the curve-of-growth, we need to assume an intrinsic 
profile of the target absorption band. It is known that the absorp- 
tion profiles of the CO2 and CO ices change depending on the 
temperature and compositions of the ices (Gerakines et al. 1999; 
Chiar et al. 1998). By considering the observed variation in ice 
profiles around massive YSOs, we choose appropriate profiles 
for the CO2 and CO ice and estimated the error of our analysis 
(see App. A for details). We adopted the band strengths of the 



where F mrm = 1 - F obs /F c , 



(3) 
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Fig. 2. Continued 

Table 4. Column density of ices and properties of YSOs 



4.0 4.1 4.2 4.3 4.4 4.4 4.5 4.6 4.7 4.8 

Wavelength [am] 



Number N(H 2 0) 


N(C0 2 ) 


N(C0 2 )/N(H 2 0) 


N(CO) 


N(CO)/N(H 2 0) 


EW(Bra-) L, of 


[10 17 cm" 2 ] 


[10 17 cm" 2 ] 


[%] 


[10 17 cm" 2 ] 


[%] 


[A] [10 3 xL o ] 



ST1 


22.5 


±0.96 


5.66: 


±0.64 


25 


3.62 ± 0.93 


16 


100 ±7 


53 ± 18 


ST2 


11.9 


± 1.0 


3.64: 


±0.38 


31 


<0.5 


<1 


343 ± 12 


182 ± 9 


ST3 


38.1 


±3.9 


10.6 


± 1.6 


28 


3.17 ±0.92 


8 


21.3 ±6.3 


9.9 ± 1.0 


ST4 


19.6 


±0.83 


7.41 : 


±0.79 


38 


2.71 ± 1.2 


14 


17.4 ± 5.1 


4.0 ± 0.4 


ST5 


26.4 


± 1.2 


10.1 : 


±0.79 


39 


2.33 ± 0.82 


9 


119 ± 7 


66 ± 12 


ST6 


45.6 


± 16° 


17.1 


± 3.0 


37 


<13 




44.3 ± 21 


6.6 ± 2.5 


ST7 


27.8 


±0.95 


19.5 


± 1.8 


70 


<4 


<14 


15.8 ±6.1 


20 ±5 


ST8 


20.1 


± 1.2 


9.34 


± 1.1 


47 


2.60 ± 0.75 


13 


64.1 ± 6.2 C 


46 ± 16 rf 


ST9 






1.84: 


± 0.17 




<2 






54 ±22 


ST10 


39.8 


± 5.0 


14.5 


± 2.1 


36 


7.72 ± 1.5 


19 


<10 


32 ± 14 


ST11 


8.9 


±0.3 


0.9 ± 


0.23 


10 


<0.5 


<1 


116 ± 6. 


338 ± 27 


ST12 


2.34 


±0.12 


0.9 ± 


0.15* 


39 






<7 


23 ± 3 d 



"Derived from low-resolution spectrum of Shimonishi et al. (2008), but re-calculated after applying our own spectral reduction pipeline 
*Derived by Oliveira et al. (2009) 

c Detection of a strong [OIII] emission line at 88 /jm is reported by van Loon et al. (2010a). 

^Estimated luminosities 95 and 62L Q from van Loon et al. (2010a), on the basis of far-infrared data, for ST8 and ST10, respectively. 
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Fig. 2. Continued 



C0 2 and CO icestobe7.6xlCT 17 and l.lxlCT 17 cm molecule -1 , 
respectively (Gerakines et al. 1995). 



Gas-phase CO2 shows narrow absorption lines around the 
wavelength range of the CO2 ice feature, and it is difficult to re- 
solve it with our spectral resolution. However, the contribution of 
the gas-phase CO2 is generally small toward YSOs (Gerakines 
et al. 1999; Nummelin et al. 2001), and we assume that it has 
a negligible effect on the derived equivalent widths, except for 
ST12, which shows a strong absorption band of the gas-phase 
CO2. On the other hand, the gas-phase CO ro-vibrational lines 
sometimes appear in the spectra of YSOs and are blended with 
the CO ice feature with our spectral resolution (Chiar et al. 
1998). Unresolved gas-phase CO absorption bands are expected 
to cause shallow and broad absorption features on either or both 
sides of the CO ice feature. Therefore, in cases where these 
bands are seen on the sides of the CO ice feature, we subtracted 
them by assuming a Gaussian profile. The derived column den- 
sities are listed in Table 4. 



4.2.3. CH3OH ice, °C0 2 ice, and XCN 

In this study, we did not estimate the column densities of these 
minor ice species because the observed strengths of these ice 
features are too weak to discuss the column densities with the 
present S/N. In addition, it is difficult to distinguish the 13 OC>2 
ice and XCN features from the unresolved CO gas absorption 
bands by our spectral resolution. Observations of these minor 
ice features with higher spectral resolution and better S/N are 
challenges for the future. 



5. Discussion 

5.1. Large C0 2 /H 2 ice ratio in the LMC 

The derived column densities of the H2O and CO2 ices are plot- 
ted in Figure 3. The average CO2/H2O ice column density ra- 
tio of our samples is calculated to be 0.36 + 0.09. The present 
samples are massive YSOs and should be compared to Galactic 
massive samples. Thus the column densities of Galactic massive 
YSOs (CO2/H2O ~ 0.17 ± 0.03, Gerakines et al. 1999; Gibb 
et al. 2004) are also plotted in the figure for comparison. It is 
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reported that Galactic quiescent molecular clouds show a simi- 
lar C0 2 ice abundance (~ 0.18 ± 0.04, Whittet et al. 2007). As 
can be seen in the figure, the typical abundance 3 of the CO2 ice 
toward the LMC's massive YSOs is larger than for Galactic mas- 
sive YSOs. Also, the scatter of the CO2/H2O ratio is larger than 
for Galactic samples. 
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Fig. 3. CO2 ice vs. H 2 ice column density in units of 10' 7 cm~ 2 . Open 
squares with error bars represent the results of this study. Filled squares 
represent those of Galactic massive YSOs (Gerakines et al. 1999; Gibb 
et al. 2004). An average CO2/H2O ice ratio of our samples (0.36) is 
plotted as a solid line. Dotted and dashed lines represent CO2/H2O ra- 
tio ~ 0.32 and 0.17, which correspond to the typical CO2/H2O ratio 
of Galactic low- and intermediate-mass YSOs, and massive YSOs, re- 
spectively. The dotted line also corresponds to the CO2/H2O ratio of the 
LMC's massive YSOs derived by Oliveira et al. (2009). 



Shimonishi et al. (2008) reported that the typical C0 2 /H 2 
ice ratio in the LMC is 0.45 + 0.17 based on the low-resolution 
spectra of the LMC's YSO (targets included in the present sam- 
ples). The present result is consistent with the CO2/H2O ratio of 
Shimonishi et al. (2008) within the uncertainties, but it is slightly 
lower than our previous result. It is mostly attributed to the effect 
of the PAH 3.3/mi emission, which was not resolved and thus 
not correctly counted in the previous study. Oliveira et al. (2009) 
have estimated the column densities of the CO2 ice for some 
common objects with the present targets based on the 15.2yum 
feature and also revised some of the column density estimates 
of Shimonishi et al. (2008). The CO2 ice column densities esti- 
mated in Oliveira et al. (2009) are slightly lower than our esti- 
mate in Table 4, but generally consistent within the uncertainties. 
They also suggest that the CO2 ice abundance of the LMC sam- 
ples is higher than in massive YSOs in our Galaxy (CO2/H2O ~ 
0.32, Fig. 3). 

The high CO2/H2O ratio in the LMC indicated by the previ- 
ous studies is now confirmed with the increased number of sam- 
ples. In addition, the uncertainties of the derived column densi- 
ties decreased from those of Shimonishi et al. (2008) thanks to 
the higher S/N, the higher spectral resolution, and the improved 
treatment of the curve-of-growth method. It is also worth men- 
tioning that a relatively high C0 2 /H 2 ratio of 0.32 + 0.02 is 
observed toward Galactic low- and intermediate- mass YSOs 



3 The "abundance" of a given ice is generally defined as the ratio of 
the column density relative to the H2O ice since it is the most abundant 
ice species. 



(Pontoppidan et al. 2008). Their samples also show a large scat- 
ter in the CO2 ice abundance as seen in the present samples. A 
similar trend in the CO2 ice abundance between the LMC's mas- 
sive YSOs and Galactic low- and intermediate-mass YSOs may 
indicate some similarities of star formation conditions between 
these different samples, but a qualitative explanation is an issue 
for future investigation. 

The derived CO/H 2 ice ratio of the LMC's YSO ranges 
from ~0.01 to 0.2 (Table 4), and this is similar to what is ob- 
served toward Galactic massive YSOs (~0.02 - 0.2, Chiar et al. 
1998; Gibb et al. 2004). Also, the distribution range of the H 2 
ice column densities of our samples is nearly comparable to 
that of the Galactic YSOs. Therefore, we can assume that there 
should be a different physical or chemical environment that se- 
lectively enhances the production of CO2 ice in the LMC. 

It is widely accepted that the CO2 ice in the cold envelope 
of a YSO is formed by grain surface reactions since gas-phase 
reactions are not able to produce the observed abundance of the 
CO2 ice, while molecules like CO are mainly formed by the gas- 
phase reaction (e.g., Millar et al. 1991). However, the formation 
mechanism of the CO2 ice still remains unclear, and it is one of 
the key topics in astrochemistry. Laboratory experiments indi- 
cate that the CO2 ice is efficiently produced by the UV photon 
irradiation to H2O-CO ice mixtures (e.g., Gerakines et al. 1996; 
Watanabe et al. 2007). It is reported by several authors that the 
LMC has an order-of-magnitude stronger UV radiation field than 
our Galaxy owing to its active star formation and its metal-poor 
environment (Israel et al. 1986; Lee et al. 2007). Thus the strong 
UV radiation field is able to help the efficient formation of the 
C0 2 ice in the LMC. 

On the other hand, an alternative mechanism of the CO2 
ice formation has been proposed to explain the detection of the 
abundant CO2 ice toward quiescent molecular clouds, where en- 
ergy sources of UV photon are not expected (Whittet et al. 1998). 
A theoretical study suggests that the CO2 ice can also be formed 
through the diffusive surface chemistry without UV irradiation 
(Ruffle & Herbst 2001). This model suggests that a sufficient 
abundance of the CO2 ice can be produced at relatively high 
dust temperatures, but this process is quite sensitive to the dust 
temperature and other assumed conditions. Several studies re- 
port that the interstellar dust temperature in the LMC is gen- 
erally higher than in our Galaxy based on far-infrared to sub- 
millimeter observations of diffuse emission (e.g., Aguirre et al. 
2003; Sakon et al. 2006). In addition, van Loon et al. (2010a,b) 
derived the dust temperatures of YSOs in the LMC and Small 
Magellanic Cloud (SMC) and suggest that the dust temperatures 
become higher in more metal-poor environments of the SMC. 
Thus the dust temperatures of the present samples may be sys- 
tematically higher than their Galactic counterparts, which are in 
more metal-rich environments. Therefore the high dust temper- 
ature in the LMC can also be a possible cause for the enhanced 
production of the CO2 ice. From the discussion, we conclude that 
the general properties of the LMC's environments, the strong 
UV radiation field, and the high dust temperature are responsi- 
ble for the observed high CO2 ice abundance. 

One of the interstellar molecular features that likely develop 
from UV photolysis is the 4.62/mi "XCN" feature, which is 
probably due to the OCN (e.g., Pendleton et al. 1999; Bernstein 
et al. 2000). This feature provides a diagnostic of the local ra- 
diation field of objects (Spoon et al. 2003). In our galaxy, the 
XCN feature is generally very weak or absent toward YSOs and 
dark clouds, but is very strong toward a few high-mass YSOs 
and Galactic center objects (Gibb et al. 2004). Thus this feature 
could suggest that the LMC's strong UV radiation field is a dom- 
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inant factor for the high CO2 ice abundance in the LMC. Some 
of our samples show a hint of this XCN feature on the shorter 
side of the CO ice feature (ST6, 7, 10). There may also be an 
additional component seen in the spectrum of ST4 and ST8, but 
the peak wavelength seems to be slightly shorter than the XCN 
feature. No clear trends are seen in the CO2/H2O ratio of the 
objects with possible detection of the XCN. However, we can- 
not separate this feature from unresolved CO gas lines because 
the spectral resolution is not high enough, and the detection of 
XCN is still tentative, so that no definite conclusion can be drawn 
from the present spectra. Future observations should provide the 
needed resolution. 

5.2. Correlation of ice abundance with YSO properties 

As described in §2, the present samples basically possess similar 
characteristics (high-mass Class I objects). Despite this, column 
densities of the H2O, CO2 and CO ices and their ratio show large 
variations. To investigate which physical condition is dominant 
in the formation and evolution of ices, we compare the chemi- 
cal conditions of ices and the properties of individual YSOs. It 
is probable that the radiation from the central star dominates the 
circumstellar environment of a massive YSO. If the formation of 
ices is controlled by the radiation from the star, there should be 
some parameters of YSOs that correlate well with the ice abun- 
dances. 

First, we compared the strength of the hydrogen emission 
lines with the CO2 and CO ice abundance. Although it is not 
seen in the low-resolution spectra of Shimonishi et al. (2008), 
the present spectra clearly detect hydrogen recombination lines 
toward the target YSOs. It is known that massive YSOs in their 
more evolved stage ionize their circumstellar gas by stellar UV 
radiation and form compact HII regions around the central star 
(van den Ancker et al. 2000). A variation in their hydrogen line 
strength thus indicates that the present targets are in their various 
evolutionary stages and possess different radiation environments 
(especially as for UV) in their circumstellar region, as discussed 
in Seale et al. (2009). We measured the equivalent width of the 
Bra line at 4.05/im, which is the strongest hydrogen recombina- 
tion line in the wavelength range of our observations. 

Figure 5.2 shows the comparison of the equivalent width of 
the Bra line with the CO2/H2O and CO/H2O ice ratio. It seems 
that the abundance of the CO2 ice does not correlate with the 
strength of the Bra line. The objects with similar CO2 abundance 
show various Bra line strengths. This result indicates that the 
UV radiation from the central star is a less dominant factor for 
the evolution of the CO2 ice. However, in the previous section, 
we mentioned that the UV irradiation can help form CO2 ice. 
Several authors have demonstrated by simple calculation that 
the interstellar radiation field can penetrate the site where the 
C0 2 ice forms (e.g., Whittet et al. 1998; Watanabe et al. 2007). 
We therefore suggest that the external (interstellar) UV radiation 
field may play an important role in CO2 ice formation. On the 
other hand, it seems that objects with strong Bra line strengths 
(ST2, ST1 1) show a small CO ice abundance (upper limit), and 
one with the weakest Bra line (ST 10) shows the highest CO ice 
abundance. Although the CO ice abundance of relatively weak 
Bra objects is rather scattered, the result indicates that the CO 
ice around YSOs is more likely to be affected by the stellar radi- 
ation than the CO2 ice. 

A similar comparison was also made for the total luminosity 
of the YSOs, which is also a good estimate of the radiation envi- 
ronment of a YSO. The total luminosity was derived by the SED 
fitting described in §2. As seen in Fig. 5.2, no clear correlation 



is seen between the CO2 ice abundance and total luminosity of 
YSOs. The trend seen in the plot of the CO ice against the total 
luminosity is somewhat similar to what is seen in the compari- 
son with the Bra strength; i.e., the overall correlation is weak, 
but the most luminous objects in the present sample (ST2 and 
ST1 1) show a small abundance of CO. 

The low CO ice abundance toward luminous and probably 
evolved (indicated by strong Bra emission) YSOs can be ex- 
plained by the low sublimation temperature of the CO ice. The 
sublimation temperature of the pure CO ice is lower than that of 
the pure H 2 and C0 2 ices (16K, 50K, and 90K, for CO, C0 2 , 
and H2O, Tielens 2005). Although the sublimation temperature 
changes depending on the chemical compositions of ices, their 
general trend does not change. The CO ice is therefore more 
easily affected by the stellar radiation owing to its low subli- 
mation temperature. Given a similar geometry and distance be- 
tween the central star and the circumstellar dust, the dust around 
more luminous YSOs is expected to be warmer since the radia- 
tion from the central star is the dominant heating source of their 
circumstellar dust. Likewise, dust around more evolved objects 
is expected to be warmer thanks to the dissipation of the cir- 
cumstellar dust. The relatively high luminosity of ST2 and ST1 1 
suggests that these objects might be clusters. However, we re- 
gard them here as YSOs of a similar kind as in other samples, 
because their infrared SEDs are well-fitted by the single YSO 
model of Robitaille et al. (2007), as well as in other samples. We 
thus suggest that the CO ice around the present luminous sam- 
ples (ST2 and ST11) may be sublimated thanks to the intense 
radiation from the central star. The CO abundance of less lumi- 
nous objects also has a scatter, which suggests that other factors 
than the stellar radiation may be important for the CO ice chem- 
istry as for these objects. 

6. Summary 

To understand the chemical conditions of ices around YSOs in 
the Large Magellanic Cloud, we performed spectroscopic ob- 
servations toward massive embedded YSOs and the candidates 
using the IRC onboard the infrared satellite AKARI. We obtained 
near-infrared spectra (R~80, 2.5-5um) of 12 objects in total. The 
ice absorption features of H 2 0, CH 3 OH, C0 2 , 13 C0 2 , CO, and 
possibly XCN were detected in the spectra. In addition, hydro- 
gen recombination lines and PAH emission bands were detected 
toward the majority of the targets. 

Compared to our previous study (Shimonishi et al. 2008), the 
increased number of LMC's YSO samples, the higher spectral 
resolution, and the improved treatment of the curve-of-growth 
method make it possible to derive column densities of ices and 
their typical abundances more accurately. Also, the detection 
of hydrogen emission lines and the SED analysis of individ- 
ual YSOs allow discussion of the correlation between ice abun- 
dances and YSO properties. 

We derived the column densities of the H2O, CO2 and CO 
ices by fitting the laboratory ice spectra (for 3.05//m H2O ice fea- 
ture) and by the curve-of-growth method (for 4.27//m CO2 and 
4.67/im CO ice features). As a result, it is shown that the typi- 
cal CO2/H2O ice ratio of our samples (~0.36 ± 0.09) is higher 
than that of Galactic massive YSOs (-0.17), while the CO/H 2 
ice ratio is comparable to Galactic ones. The systematical dif- 
ference in the CO2 ice abundance around the LMC's massive 
YSOs, which was suggested by previous studies (Shimonishi 
et al. 2008; Oliveira et al. 2009), is now confirmed with the 
new near-infrared data. It is also shown that the scatter of the 
CO2/H2O ratio in the LMC is larger than for Galactic ones. 



10 



Takashi Shimonishi et al.: Ices around Embedded YSOs in the LMC 




100 10 
Equivalent Width (Bra) [A] 



0.8 
£"0.6|- 

z 

0.4 h 

/— > 

CM 

8 0.2 h 



0.0 



—I 1 — I- 



* ■ 




10 



100 10 
Total Luminosity [L Sun xio 3 ] 



100 



Fig. 4. The equivalent width of the hydrogen recombination line (Bra) vs. CO2/H2O ratio (upper left) and CO/H2O ratio (upper right). And a total 
luminosity of sample YSOs vs. CO2/H2O ratio (bottom left) and CO/H 2 ratio (bottom right). A YSO number given in Table 1 is plotted for the 
objects discussed in the text. A value of ST6 is not plotted due to its large uncertainties. ). 



Although the formation mechanisms of the CO2 ice is still 
unclear, we suggest that the strong ultraviolet radiation field 
and/or the high dust temperature in the LMC are responsible for 
the observed high abundance of the CO2 ice. Our result provides 
evidence that YSOs in the metal-poor LMC maintain different 
chemical conditions in their circumstellar environments. 

A relatively high CO2 ice abundance (~0.32) and a large 
scatter are also reported toward Galactic low- and intermediate- 
mass YSOs. A similar trend of the CO2 ice abundance between 
the LMC's massive YSOs and Galactic low- and intermediate- 
mass YSOs may indicate some commonality in these objects, 
but a qualitative explanation considering both the structure of a 
YSO envelope and the ice chemistry is required in the future. 

The present YSO samples possess similar characteristics in 
their infrared SEDs, but they show a large variation in the CO2 
and CO ice abundances and in their near-infrared emission fea- 
tures. We compared the CO2 and CO ice abundances with the 
strength of the hydrogen recombination line and total luminosity 
of the sample YSOs. It is shown that the CO2 ice abundance and 
these YSO characteristics do not show a correlation. The present 
result indicates that the internal stellar radiation does not play an 
important role in the evolution of the CO2 ice around a mas- 
sive YSO. We thus need further investigations to shed light on 
the dominant physical factors in the CO2 ice formation. Similar 
to the CO2 ice, a clear correlation is not seen between the CO 
ice abundance and the YSO properties. However, it is suggested 
that the CO ice abundance of luminous and probably evolved 
samples are significantly lower than other samples. We infer that 



more volatile molecules like CO are susceptible to the effect of 
the stellar radiation. 

The present results revealed the presence of a massive YSO 
that shows a CO2/H2O ratio of up to ~0.7. If there is also a 
chemical difference in the LMC's lower mass YSOs, it will in- 
directly tell us the chemical diversity of subsequently-formed 
planets since the interstellar ices are thought to be an origin of 
planetary system ices, as discussed in Ehrenfreund & Schutte 
(2000). The future infrared telescope missions SPICA, JWST, 
and SOFIA will allow better understanding of the diversity of 
solid materials in the universe. 
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Appendix A: The curve-of-growth method 

We calculated column densities of unresolved ice features by us- 
ing the equivalent width. This calculation technique is called the 
"curve-of-growth" method, and it is traditionally used to study 
gaseous absorption lines (Shore 2002). Since the optical depth is 
the logarithm conversion of the normalized flux (1 - F /, s /F e ), the 
peak optical depth derived from unresolved features is not a true 



Table A.l. Ice profiles used for calculating the curve-of-growth 



Molecule 



Ice Profile 



C0 2 
CO 



H 2 : C0 2 ice mixture at 14K 
Gaussian profile (FWHM=6.89 cirT 1 , A cenler =4.67fim) 



value. Thus, in case of our low-resolution spectra, we will sig- 
nificantly underestimate a column density if it is derived from 
the simple integration of the optical depth as expressed in Eq. 
(1), and this is why we need to use the curve-of-growth method. 
This method was applied for the first time to absorption features 
of ices in Shimonishi et al. (2008), but we overestimate the un- 
certainties and a detailed description of the method was not given 
because of space limitations. In this appendix, we describe our 
improved treatment of the curve-of-growth method. 

First, we convert the measured equivalent width into a col- 
umn density. To do this, we need to assume intrinsic absorption 
profiles of target ice features. Here, they are the 4.27//m CO2 ice 
and 4.67/im CO ice feature. For CO2, we considered the typi- 
cal compositions and temperature of the interstellar CO2 ice and 
chose an appropriate laboratory spectrum as shown in Table A. 1 
(Ehrenfreund et al. 1996, data taken from the Leiden Molecular 
Astrophysics database). For CO, it is reported that several com- 
binations of laboratory ice profiles can account for the same ob- 
served spectrum and the fitting degenerates (Pontoppidan et al. 
2003). We therefore used a single Gaussian profile that is close 
to the typically observed CO ice feature around massive YSOs 
(Chiar et al. 1998). The assumed profiles of the C0 2 and CO 
ices features are summarized in Table A.l. By using these pro- 
files, we calculated both equivalent widths and column densities 
for the various optical depths. The plot of the derived equiva- 
lent widths vs. column densities is shown in Fig. A.l, and this is 
called the curve-of-growth. The equivalent widths can be mea- 
sured from observations, and these plots are used to derive col- 
umn densities of ices; however, the uncertainty becomes larger 
for the higher column densities owing to the saturation effect of 
the curve-of-growth. 

It is known that the actual profile of the observed ices in star- 
forming regions is not uniform since they vary depending on 
the compositions and temperature of ices (Gerakines et al. 1999; 
Gibb et al. 2004). The derived column densities depend on the 
assumed profiles, so we need to check the validity of our method 
and to evaluate the uncertainties caused by the method. 

We first prepared 2-5/j.m ISO-SWS spectra of embedded 
massive YSOs 4 (Sloan et al. 2003), whose ice column densities 
are already presented in Gibb et al. (2004). The spectral resolu- 
tion of these data is high enough to resolve narrow ice features 
(R ~500-1000). Then their spectral resolution was degraded to 
the resolution of the AKARI IRC NG mode (R -80) by convolv- 
ing the point spread function of the IRC. Next, we derived col- 
umn densities of the CO2 and CO ice by these converted spectra 
using the same method as that described in §4.2.2. The derived 
column densities and the results of Gibb et al. (2004) are plotted 
in Fig. A. 2. 

The figure indicates that the column densities derived from 
our curve-of-growth method are basically in good agreement 
with the results of Gibb et al. (2004), which derived the column 
densities from the high-resolution spectra. But a systematic dif- 
ference and scatter are also seen in the figure. These are mainly 



Data taken from http://isc.astro.cornell.edu/ sloan/library/swsatlas/atlas.html 
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N(C02) by profile fitting N(CO) by profile fitting 

Fig. A.2. Comparison of the CO2 and CO ice column densities derived from the low-resolution spectra by the curve-of-growth and those derived 
from the high-resolution /50-SWS spectra by the laboratory ice fitting (Gerakines et al. 1995; Gibb et al. 2004). The axes are expressed in units 
of 10 17 cnr 2 . Solid and dashed lines indicate the result of the linear fit and its one sigma error. The labeled numbers in the figure correspond to 
the name of the objects that are used in this analysis as follows; 1. Mon R2 IRS3, 2. Orion IRC2, 3. S140, 4. Elias 29, 5. AFGL 2136, 6. W33A, 
7. NGC7538 IRS9. 



caused by the difference between the observed and assumed ice 
profiles. The presence of weak emission or absorption features 
may also be responsible for these differences. We obtained the 
following equations from a straight-line fit to the data points in 
Fig. A.2. 

N(C0 2 ) = N C 04 x (0.98 ± 0.038) (A.l) 
and 

N(CO)=Ncq4X (1.08 + 0.065), (A.2) 

where Nco4 represents the column density of each ice presented 
in Gibb et al. (2004). In Shimonishi et al. (2008), we adopted 
20% uncertainties to account for the above systematic differ- 
ence in the derived column density of the CO2 ice and added 
to the observational uncertainties. The uncertainties were, how- 
ever, overestimated as expressed in Eq. (A.l). In this study, the 
column densities derived from the curve-of-growth method are 
divided by the slope of the corresponding line to correct the sys- 
tematic difference in the derived column densities. In addition, 



a 1-cr fitting error of the slope is considered as an uncertainty 
of our curve-of-growth method and added to the uncertainties in 
the observation and the fitting. 



